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Abstract we examine the effect of thermal conduction 
on the observational properties of a super critical hot 
magnetized flow. We obtained self-similar solution of 
a magnetized disc when the thermal conduction plays 
an important role. Follow of our first paper (Ghasem- 
nezhad et al. 2012 (hereafter GKA12)) we have ex- 
tended our solution on the observational appearance of 
the disc to show how physical condition such as thermal 
conduction, viscosity, and advection will change the ob- 
served luminosity of the disc, Continuous spectra and 
surface temperature of such discs was plotted. We ap- 
ply the present model to black-hole X-ray binary LMC 
X-3 and narrow-line seyfert 1 galaxies, which are sup- 
posed to be under critical accretion rate. Our results 
show clearly that the surface temperature is strongly 
depends on the thermal conduction, the magnetic field 
and advection parameter. However we see that thermal 
conduction acts to oppose the temperature gradient as 
we expect and observed luminosity of the disc will re- 
duce when thermal conduction is high. We have shown 
that in this model the spectra of critical accretion flows 
strongly depends on the inclination angle. 
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1 INTRODUCTION 

Accretion discs are found around several kind of as- 
trophysical objects: from young stars (protoplanetary 
discs) to supermassive black holes in AGNs. The dy- 
namics of these discs is however poorly understood. Ac- 
cretion process is penetrated itself in various ways de- 
pending on that how angular momentum of the discs it 
carries away and how the its energy is dissipated. Such 
systems typically under the effect of several physical 
process. It is known that, on average, matter moves 
inward resulting in the accretion of material onto the 
central object and powers the most luminous objects 
in the universe by converting gravitational potential 
energy into highly energetic radiation. Accretion disc 
forms when gaseous matter, usually an assembly of free 
electrons and various types of ions, spirals onto a cen- 
tral gravitating body by gradually losing its initial an- 
gular momentum effectively in a timcscales compatible 
with observations, as a result of viscous and magnetic 
stresses. The energy released by accretion discs de- 
pends on the mass accretion rate as well as spin and 
mass of central black hole. There is a maximum value 
for luminosity (Eddington luminosity) of accretion disc 
which gravity is able to exceed the outward pressure of 
radiation. 

The theory of black hole accretion disc has been de- 
veloped rapidly since 1970s and brings a lot of successes 
in describing high-energy astronomical phenomena (See 
for a review Kato et al. 2008). As for thermally stable 
black hole accretion disc models, there are only three 
classes have been developed. The fist one is the stan- 
dard Shakura & Sunyaev disc, which is optically thick 
and geometrically thin, Keplcrian rotation with very 
small infall velocity, and truncated at the marginally 
stable orbit. The standard disc was studied for the 
Newtonian case (e.g., Shakura & Sunyaev 1973) and for 
the relativistic case (e.g., Novikov & Thorne 1973). The 



second one is the optically thin ADAFs (advection dom- 
inated accretion flows) that works very well for small 
mass-accretion rates (Narayan & Yi 1994) while for su- 
percritical accretion discs (M ~ M crit ) slim discs was 
introduces (Abramowicz ct al. 1988). M crit is the criti- 
cal mass accretion rate, defined by M C Ht = — ?• Le and 
c are the Eddington luminosity and light speed, respec- 
tively. Therefore this disc brightness about Edding- 
ton or supcr-Eddington luminosity and the geometrical 
thickness cannot ignore (— < 1) (Abramowicz et al. 
2010). The Slim disc models describes supercritical disc 
accretion flow, so they are identical to each other. Re- 
cently the slim disc model began to attract considerable 
attention, because of its success to reproduce spectra 
of accreting binaries, the spectral behaviour of several 
black hole binaries cannot be explained by the tradi- 
tional standard disc, but is well reproduced by slim disc 
models (Kubota 2001). LMC X-3 is black hole (BH) bi- 
nary system in the large magellanic cloud (LMC) at a 
distance of 48.1 kpc (derived from Orosz et al. 2009). 
Cowley et al. (1983) stablished LMC X-3 as a BH con- 
didate with a BH mass 7M© < M B h < 14M . LMC 
X-l has a soft spectra and a low absorption column 
density along the line of sight. These properties make 
LMC X-3 an ideal sample for testing our understand- 
ing of black hole accretion disc physics. Slim disk model 
describes accretion flows at high luminosities, while re- 
ducing to the standard thin disk in the low luminusity 
limit. The most important feature of super critical ac- 
cretion flows is the presence of photon trapping as the 
advection of radiation entropy. The Slim disc will ex- 
ceed Eddington limit, because the accretion rate can 
be larger than Eddington accretion rate (Mineshige et 
al. 2000, Ohsuga et al. 2002). The basic stucture and 
spectral properties of the slim disks have been exam- 
ined by several researchers: Watarai & Fukue (1999) 
showed the existence of radiative winds from the su- 
perdisk, Mineshige et al (2000) applied slim disks to 
the narrow-line Seyfert 1 galaxies. Fukue (2000) ex- 
amined the basic properties of super critical accretion 
discs and derived the relations between the observables 
and the model parameters. 

Several physical process should be take into account 
when we are studying black hole accretion disc. Mag- 
netic field palys a great role on the structure of the 
disc since in hot accreting flow, the temperature is so 
high that the accreting materials are ionized. Many 
researcher have been attack the problem to solve the 
magnctohydrodynamics (MHD) equations of hot mag- 
netized flow analytically (Akizuki & Fukue 2006, Ab- 
bassi ct al. 2008, Ghanbari et al 2007, Shadmehri 
& khajenabi 2005). Kaburaki (2000) has presented a 
set of analytical solutions for a hot accretion flow in a 



global magnetic field. Ghanbari et al. (2007) have pre- 
sented a set of self-similar solutions for two-dimensional 
(2D) viscous-resistive ADAFs in the presence of a dipo- 
lar magnetic field of the central accretor. They have 
shown that the presence of a magnetic field and its as- 
sociated resistivity can considerably change the picture 
with regard to accretion flows. 

One of the largely neglected physical process in the 
physics of accretion discs is thermal conduction; Hot 
accretion flows have high temperature, so the internal 
energy per particle is high. This is one of the rea- 
sons why advection cooling overcomes radiative cool- 
ing. For the same reason turbulence heat transport 
by thermal conduction is non-negligible in heat bal- 
ance in the disc (Kato ct al. 2008). while some re- 
cent observations (Loewenstein et al. 2001; Di Matteo 
et al. 2003; Ho et al. 2003) of the hot accretion flow 
around active galactic nuclei indicated that it should 
be on collision-less regime. So thermal conduction has 
a great role in energy transport in the accreting mate- 
rials in a hot accretion disc where they arc completely 
ionized. The wcckly-collisional nature of hot acccrction 
flows has been noted before (Mahadevan & Quataret 
1997). Since thermal conduction act to oppose the for- 
mation of the temperature gradient that causes it, that 
might expect that the temperature profile in a thermal 
conducting disc should be differ from the case of the 
disc was not under the influence of thermal conduction. 
Shadmehri (2008), Abbassi et al. (2008, 2010), Tanaka 
& Menou (2006) have studied the effect of hot accre- 
tion flow with thermal conduction in a semi-analytical 
method; physics of such systems have been studied in 
simulation models (e.g. Sharma et al. 2008; Wu et 
al. 2010). Abbassi ct al. (2008) have shown that for 
this problem there are two types of solutions ; high 
and low accretion rate. Ghasemnezhad et al. (2012) 
have shown, by putting an extra physical condition 
have shown that the high accretion rate solutions of 
Abbassi et al. (2008) are not exactly correct as long as 
some of the low accretion rate solutions arc not phys- 
ical meaning. This extra condition make a acceptable 
parameter space which in this space the self-similar so- 
lutions has reasonable physical behaviour. Following 
GKA12 we are investigating the observational appear- 
ance a hot magnetized flow under effect of thermal con- 
duction using Fukue (2004) framework in a parame- 
ter space presented by GKA12. We developed Fukue 
(2004) solutions for narrow-line seyfert 1 galaxies by 
adding the thermal conduction and toroidal magnetic 
field effect. We suppose a moderately massive central 
black holes ~ 10 5_6 M Q under critical accretion rate 
IMoyr- 1 (e.g.,Minishige et al. 2000). 

This paper is organized as follows: Section 2, we 
present the equations of magnetohydrodynamics as the 



basic equations and assumptions. Self-similar solutions 
are presented in section 3. In section 4, the continuum 
spectra are calculated and we regard general properties 
of supercritical discs and ultimately we show the results 
in section 5. 



2 The Basic Equations 

For constructing the model, we used a set of coupled 
differential equations describing the law of conserva- 
tion for a steady state, axi-symmetric (^ = -J^ = 0) 
supercritical accretion disc. In cylindrical coordinates 
(r, if, z), we used vertically integrate the flow equations, 
also, we suppose that all variables are only a function 
of r. We ignore the relativistic effect and we use New- 
tonian gravity. We adopt a-prescription for viscosity of 
accreting flow. The magnetic field was considered with 
toroidal configurations. 

The MHD equations are as the same as (Akizuki & 
fukue 2006, Abbassi et al. 2008, GKA12): 
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parameter in our model. Following (GKA12) radiative 
pressure is ignored but it should take into account in 
the case of slim disc to have more realistic picture. As 
the same as Akizuki & Fukue 2006, Abbassi et al 2008 
and GKA12 we can chose two cases: Case 1: when 
the pressure is assumed to be the gas pressure (thermal 
pressure). Case 2: when the pressure is assumed to be 
the magnetic pressure plus the gas pressure. So in the 
case 2 we replace a whit a(l + B) in all of the equations. 
For simplicity in this investigation we use the case 1. 

We will follow how will change is not here the dynam- 
ical and observational properties of the disc will change 
when the /? varies is a physically reasonable interval. 
In the energy equation we have cooling and heating 
term in the disc. We assume the generated energy due 
to viscous dissipation plus the heat conducted into the 
volume are balanced by the advection cooling. The 
second term on right hand side of the energy equation 
represents energy transfer due to the thermal conduc- 
tion where F s = 5$ s/ oCg is the saturated conduction 
flux (Cowie & Makee 1977). Dimensionless coefficient 
<fr s is less than unity. 

The induction equation (the final one) represents the 
interaction of fluid and magnetic field. Where B v is the 
field scaping/creating rate due to magnetic instability 
or dynamo effect. 
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where V r is the accretion velocity (V r < 0) and E = 
2pH is the surface density at a cylindrical radius r, p 
is the mass loss rate per unit volume, H is the disc 
half-thickness. The azimuthal equation of motion was 
integreated over z. a is the viscous parameter, J7(= 
— ) and ilk are the angular velocity and the Keplcrian 
angular speed respectively. 

In this equation we have introduced 8 = rT ag = 

^gas 

i(^p) 2 which is magnetic field pressure over the gas 
pressure which indicates the importance of magnetic 
field in the dynamics of accretion flow and it is a free 



Fig. 1 The surface temperature of LMC X-3 as a func- 
tion of dimensionless radius ( — ) for several values of (top 
panehmagnetic field (/3) , middle panebadvetion parameter 
(/), bottom panebthermal conduction {4>s))- 



3 Self-Similar Solutions 

To solve the MHD equations which they are introduced 
in last section and have physical interpretation of a hot 
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Fig. 2 Continuum spectra of LMC X-3 for several values of 
(top panebthermal conduction (<f> s ) , middle panehadvction 
parameter (/), bottom panebmagnetic field (/3)). 



accretion flow under effect of thermal conduction and 
toroidal magnetic field, we seek self-similar solutions 
of the above MHD equations. The self-similar method 
is familiar from its wide applications in astrophysics 
and it's abilities to solve the full set of MHD equa- 
tions. Usefulness of the self-similar method is simplify- 
ing multi-dimensional problems of geometric and kine- 
matic complexity, making them analytically tractable. 
The self-similar method is not able to reproduce the 
exact solution of accretion flows, because no boundary 
conditions have been taken into account. However, as 
long as we are not interested in the behaviour of the 
flow near the boundaries, such solutions are very use- 
ful. Following GKA12 solutions we have introduce the 
velocities as follows 
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Fig. 3 Observed flux of LMC X-3 with an inclination an- 
gles, i = 66 for cj>s = 0.004 (solid lines), cf> 3 — 0.006 (dashed 
lines) and </>„ — 0.008 (dotted lines) 
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Fig. 4 Observed flux of LMC X-3 with two inclination 
angles (i = 0,i = 66), <j> 3 = 0.004 (solid lines), 4> s = 0.008 
(dashed lines) 



and constant ci, C2 and C3 are determined later from 
the magnctohydrodynamic equations. From Eq. 4 we 
will obtain the disc half-thickness H as: 



— = Jc 3 (l + /3) = tan a 
r 



(11) 



If we assume a power law form for the surface density 
E as: 



E = E r s 

we will obtain: 

p = par 2 
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Substituting the above self-similar transformation in 
to the MHD equations of the system, we '11 obtain 
the following system of dimensionlcss equations, which 
should be solve to having c±, C2 and c 3 : In case 1: 
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which is determine the behaviour of the radial flow is 
specified by the input parameter of the fluid such as a, 
$ s , B and /. If we c\ achieve from the above equation 
then we can easily have ci and ci by: 
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After algebraic manipulations, we obtain a forth or- 
der algebraic equation for c\\ 

D 2 cj + 2DBcl + (B 2 -2D)c 2 1 -(2B + A 2 )c 1 + l = (20) 
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which it is 



= — i we have p'o 
correspond to no mass loss or wind in the hot magne- 
tized flow, Abbassi et al. (2008). In this work we focus 
on constant accretion rate ( no outflow ) case but for 
s > — 4 there is wind or outflow and accretion rate is 
not constant. The observation evidence shows that the 
outflow can exist in ADAF (like sgrA*) and Slim disc ( 
like Scyfert 1 sources). The properties of outflow exam- 
ined by Xie & Yuan (2008), Yuan, Wu, Bu 2012, Yuan, 
Bu, Wu 2012 for ADAF, and Ohsuga, Mineshige, et al. 
2009 for ADAF, SSD and Slim discs . The outflow will 
affect the dynamics, and the radiative output of the 
accretion disc ( Xie & Yuan 2012). It is clear that C\ 



c 3 = ci(- 
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(25) 
3(s + 1) 

GKA12 by adding one more constrain on solution 
presented by Abbassi et al 2008, have shown the this 
new physical condition places a constraint on the phys- 
ically valid parametric space. Their calculations arc 
carried out in a range of viscosities, thermal conduc- 
tion, and magnetic field parameters which they have 
acceptable physical behaviour. GKA12 have presented 
dynamical behaviour of such flows by plotting ci, C2, C3. 
In this paper we want to further explore this system to 
investigate the possible role of thermal conduction and 
magnetic field on the observational appearance of hot 
magnetized flow. In the next section we will follow this 
idea. 



4 radiation properties 

We have a complete set of equations which describe the 
global dynamics of hot magnetized flow. Using the re- 
sults obtained in the previous section we will able repro- 
duce observational appearance of the hot magnetized 
flow. The surface flux and disc luminosity of the discs 
are derived as follow: By assuming a dominance of the 
radiation pressure, wc can write the height-integrated 
pressure II(= Sc s 2) and the averaged flux F as: 
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where a is the Stefan-Boltzman constant. The opti- 
cal thickness of the disc as the same as height-integrated 
pressure and averaged flux has a radial dependency as 



r = — feS = — kYiQr s 
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where k(= 0Acm 2 q 1 ) and r = IMS AxT are the 

electron-scattering opacity and Schwarzschild respec- 
tively. 

Hence, the effective flux and the effective tempera- 
ture of the disc surface become: 
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where Le = Anc^j^- is the Eddington luminosity. 

By integrating these equations radially we will have 
the disc luminosity as: 
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2 x 10 5 r s (Fukue, 2000). This 



equation represents, the disc's luminosity which is af- 
fected by magnetic field explicitly through the (3, but 
it would be affected by thermal conduction, advection 
parameter and viscosity through the C3 implicitly. 

The simplest spectral model we can construct for an 
accretion disc if we assume the inflowing matter in a 
thermal equilibrium. We have assumed the accretion 
flow radiates away locally like a black-body radiation. 
We expect that the hottest regions at small r to radi- 
ate at the highest energies. This simple spectral model 
is referred to as the multi-colour disc. For the present 
purpose, we assume that the disc surface radiates black- 
body radiation B v with temperature T e ff(r). Then we 
can equate the disk radiative flux F(r) to a black body 
flux aT*ff{r) where T e ff(r) is the surface temperature. 
Since aT^^ir) = J TrB v (T e ff(r))dv, We can equate 
radiative flux F v (r) = nB v {T e ff{r)). Then, the con- 
tinuum spectrum (luminosity per frequency) L„ can be 
calculated (see , e.g., Kato et al. 1998) by: 
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where a factor 2 comes from that the radiation will ra- 
diated from both side of the disc. We will integereated 



luminocity in a resonable interval: (r c 
We have used black-body function: 
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The disc is believed to be fundamentally antisym- 
metric. When viewed at an angle other than face-on 
(from above or below the plane exactly), the circular 
disc appears elliptical. The greater the ellipticity, the 
more the galaxy's disc in inclined with respect to our 
line of sight. Here we are studding of disc viewed at 
different inclinations. Observed flux from an accretion 
disc depends on the distance D{— A8.1kpc) and it's in- 
clination angle i as: 
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This equation can be used for standard discs 
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Fig. 5 The surface temperature of narrow-line seyfert 
1 galaxies as a function of dimensionless radius (— ) for 
several values of (top panekmagnetic field (/3) , middle 
panehadvetion parameter (/), bottom panekthermal con- 
duction (0 S )). 



5 Results 

In this section we apply the present model of hot mag- 
netized accretion disc to two astronomical objects, in- 
cluding an accretion disc around a black hole, which 
are supposed to be under critical accretion: LMC X-3 
and narrow-line seyfert 1 galaxies. 

LMC X-3 is a powerful source of X-rays located in 
the Large Magellanic Cloud (LMC). This X-ray source 
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Fig. 6 Continuum spectra of narrow-line seyfert 1 
galaxies for several values of (top panel:thermal conduc- 
tion ((f> a ) , middle panebadvction parameter (/), bottom 
panehmagnetic field (/?)). 



is associated with a binary system with an orbital pe- 
riod of 1.7 days (close binary). The visible component is 
a main sequence B3 star whose shape has been severely 
distorted by the gravitational field of its companion. 
Although not unambiguous, the mass of the compact 
object is estimated to be order of 10 solar masses and 
more likely is considerably higher, making this one of 
the best black hole candidates. We used LMC X-3 val- 
ues to test our model. 

In figure 1 the surface temperature (T e //) is plot- 
ted as a function of the dimcnsionlcss radius ( — ) for 
LMC X-3. It is obvious that the surface temperature 
decreases as — increases. The top panel shows the ef- 
feet of magnetic field on the surface temperature. We 
see that the surface temperature decreases by adding 
magnetic field parameter f3. In the middle panel we 
show that surface temperature increases by increasing 
advection parameter. And finally in the lower panel 
we have plotted surface temprature as a function of ra- 
dial distance for several values of thermal conduction 
parameter, <p s . Slope of the temperature gradient de- 
creases when we add the role of thermal conduction, 
by adding <j) s , as expected. This is acceptable because 
the increasing of thermal conduction causes that the 
more energy transfers to outer of disc and so reduce 
the temperature gradient. To having this panel we just 
used physically allowed parameter space which was in- 
troduced by GKA12. 

The radiation spectrum of the hot magnetized flow 
is shown in figure 2 for several values of thermal con- 
duction, advection parameter and magnetic field from 
top to bottom. The central black-hole mass is fixed 
as 10 Mq, while the mass accretion-rate M is equal 
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Fig. 7 Continuum spectra of narrow-line seyfert 1 galaxies 
with m = lO 6 m0 and LMC X-3 with m = lOra© for the two 
thermal conduction parameter (f> s = 0.004 , <f> a = 0.008 



O.GMEdd- As can be seen in figure 2, the maximum of 
vL v is always of the order of the Eddington luminosity 
(L E = 10 39 ergs _1 ). This agrees with Fukue (2004) and 
observed value of LMC X-3. As we sec by increasing of 
thermal conduction, the surface temperature decreases 
and therefore luminosity of disc also decreases. But ad- 
vection parameter acts oppose to thermal conduction. 
We expect that by decreasing the advection parameter, 
/, more energy radiated away so the total luminosity 
will be higher while we see the opposite result in figure 
2 and 6, seems this behavior comes from the lack of ra- 
diation pressure. So this results indicate that in ordctr 
to have a realistic picture we should add the radiation 
pressure in the dynamics of hot magnetized flow. 

We plotted the observed flux of LMC X-3 for in- 
clination angle i = 66° for different values of thermal 
conduction in figure 3. In figure 3,4 we evaluate the 
effect of thermal conduction on the observed flux. Fig- 
ure 4 we have shown that by increasing the inclination 
angle i for two thermal conduction parameters, the ra- 
diation spectrum of the disc uF v will increase. In both 
of two inclination angles we can easily see by increasing 
of thermal conduction (f> s , the observed flux decreases. 
In this panel we have the same result with figure 2. 
The classical soft-state spectrum of the X-ray binary 
LMC X-3 in dominated by a ~ O.Bkev. This result 
agrees with Straub ct al. (2011). Davis ct al. (2006), 
have plotted the spectrum of LMC X-3 for i = 60°, 
D = 52kpc, and they have shown that the observed 
flux is approximately 1 k ™ and the peak of the spec- 
trum is located in E = lkev. 

The same panels are shown for narrow-line seyfert 
1 galaxies in figure 5 ,6. Seyfert galaxies are charac- 
terized by extremely bright nuclei, and spectra which 



have very bright emission lines of hydrogen, helium, ni- 
trogen, and oxygen. These X-ray emission may come 
from the surface of the accretion disc itself. The central 
black hole mass is 1O 6 M and M ~ !0 3 M Edd . In fig- 
ure 5 as can be seen for this object temperature is lower 
than LMC X-3. This result is according to the result 
of Abramowicz et al.(1988). We show in figure 6 the 
maximum of vL v is always of the order of the Edding- 
ton luminosity (Le = 10 44 ergs -1 ) and disc spectrum is 
peaking up v — 0.5 x 10 16 Hz. These results arc agree 
with Fukue (2004). 

In figure 7 we compared the solution for these two 
objects. We have shown that for the same thermal con- 
duction parameter when the mass of BH decreases, the 
spectrum is setting in high frequency and high energy 
while the luminosity decreases and we have hard-low 
state. 



steady. The real hot magnetized flow is a much richer 
medium, providing numerous examples of important as- 
trophysical phenomena on all scales. Thus, if we assume 
a two-temperature structure for the ions and electrons 
in the discs, then we may expect that the ions and elec- 
tron temperatures decouple in the inner regions where 
is so hot, which will modify the role of conduction con- 
sequently will modify the dynamical and observational 
appearance of the disc. One of the other limitation of 
this solution is the anisotropic character of conduction 
in the presence of a magnetic field which is indicated 
by Balbus (2001). 



6 SUMMERY AND CONCLUSION 



We have presented a simple study showing that self- 
similar solution of a hot magnetized flow under the ef- 
fect of thermal conduction. Purely toroidal magnetic 
field configuration is assumed. It was assumed the disc 
is axially symmetric and static with the a-prescription 
of viscosity. The rclativistic effects and self-gravity of 
the disc are ignored. The weak collisional nature of a 
hot accretion flow was confirmed by (Tanaka & Menou 
2006; Abbassi et al. 2008), so we have adopted a satu- 
rated form of thermal conduction as a possible mecha- 
nism for energy transportation. Observational appear- 
ance of such disc was investigated with a physically lim- 
ited parameter space presented by GKA12. 

GKA12 imposed condition puts some physical con- 
straint on self-similar solutions presented by Abbassi 
et al 2008. As a result, in order to have a physically 
valid solution, only a small part of the parameter space 
should be considered. After assuring the validity of our 
solutions by using the proper values of input parame- 
ters a, /3, 4> s and f we have investigated the influence 
of thermal conduction on the observational appearance 
of hot magnetized accretion flows. 

It is shown that the temperature gradient will de- 
crease when thermal conduction plays an important 
role, increases when thermal conduction plays an im- 
portant role. Radiation spectrum of the hot magnetized 
flow was plotted for two real cases: narrow-line seyfert 
1 galaxies and X-ray binary LMC X-3. The theoretical 
solution is compatible with observed x-ray spectra of 
these two cases. 

So far we have considered only the most basic of 
hot magnetized flow: one-fluid, axially symmetric, and 
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